

NASA Technical Memorandum 87307 


I NASA-TM-87307 19860022448 


Influence of Third-Degree Geometric 
Nonlinearities on the Vibration 
and Stability of Pretwisted, 

Preconed, Rotating Blades 


K.B. Subrahmanyam and K.R.V. Kaza 
Lewis Research Center 
Cleveland, Ohio 



m 23 a*6 

LANGLEY RESEARCH CENTER 
) IERA3Y, 


May 1986 


NASA 



NF01521 



E- 2988 


3 1176 01365 3762 


INFLUENCE. OF THIRD- DEGREE GEOMETRIC NONLINEARITIES ON THE VIBRATION 


AND STABILITY OF PRETWISTED, PRECONED, ROTATING BLADES 


K.B. Subrahmanyan* and K.R.V. Kaza 
National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

This report Is concerned with the study of the Influence of third degree 
geometric nonlinear terms on the vibration and stability characteristics of 
rotating, pretwisted, and preconed blades. The governing coupled flapwlse 
bending, edgewise bending, and torsional equations are derived Including third- 
deg re e geom e tric nonlinear elastic t erms by making use of the geometric non- 
linear theory of elasticity In which the elongations and shears are negligible 
compared to unity. These equations are specialized for blades of doubly sym- 
metric cross section with linear variation of pretwist over the blade length. 
The nonlinear steady state equations and the linearized perturbation equations 
are solved by using the Galerkln method, and by utilizing the nonrotating 
normal modes for the shape functions. Parametric results obtained for various 
cases of rotating blades from the present theoretical formulation are comparea 
to those produced from the finite element code MSC/NAS1 RAN , and also to those 
produced from an In-house experimental test rig. It Is shown that the spurious 
Instabilities, observed for thin, rotating blades when second degree geometric 
nonllnearltles are used, can be eliminated by Including the third-degree elas- 
tic nonlinear terms. Furthermore, Inclusion of third degree terms Improves the 
correlation between the theory and experiment. 


1NTR0DUC1 ION 

Considerable work has been done In the area of helicopter rotor blade 
dynamics to date. Several Investigators have derived the governing equations 
of motion for the helicopter rotor blades Incorporating various degrees of com 
plexity (refs. 1 to 3). It Is now established beyond any doubt that geometric 
nonllnearltles must be Included In the analysis for a fair prediction of rota- 
ting blade frequencies and stability boundaries. However, there remain certain 
questions concerning the degree to which the geometric nonllnearitles should 
be retained, and concerning the Initial assumptions In prescribing an ordering 
scheme ( ref s . 4 and 5) . 


*0n leave from NBKR Institute of Science and Technology, Mechanical Engi- 
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Associate, University of Toledo, Mechanical Engineering Department, Toledo, 
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Another area, somewhat similar to the helicopter rotor blade dynamics but 
further complicated due to the geometry, Is the advanced turboprop blade dynam- 
ics (refs. 6 and 7). Unlike helicopter rotor blades, the turboprop blades are 
more shell- like, possess variable sweep (a component of which may be considered 
as blade precone and typically of the order of 50° at blade tip for the turbo- 
prop blade) along their span, and are subjected to considerable rotational 
speeds which can cause large steady state deformations. Because of the complex 
geometry of the advanced turboprop blades, finite element methods are well 
suited for their vibration analysis. However, when finite element methods are 
used. It Is not practical to conduct parametric studies to obtain a physical 
Insight Into the various complicating effects. 

An Investigation was therefore undertaken by the present authors to assess 
the Individual and combined Influence of such complicating effects as sweep, 
Coriolis forces, pretwist, geometric nonllnearlties, and rotation on the cou- 
pled frequencies for blade configurations that are representative of advanced 
turboprop blades. The blades were modelled as beams having constant, but 
large, precone and linear pretwist. A preliminary study of this beam model 
for torslonally rigid blades (refs. 8 and 9) showed that both the linear and 
nonlinear Coriolis forces must be retained In analyzing thick blades while the 
Coriolis effects can safely be Ignored In analyzing thin blades. For a fair 
prediction of stability boundaries however, Inclusion of the torsional degree 
of freedom was found necessary. Thus, the coupled flapwlse bending, edgewise 
bending, torsion, and extension equations were addressed In reference 10. It 
was shown that for small precones and moderate rotational speeds, the frequen- 
cies produced by the beam theory with geometric nonllnearlties up through 
second degree are In excellent agreement with those produced by MSC/NAS1RAN. 
However, for large precones and large rotational speeds, the beam theory with 
second degree geometric nonllnearlties was found to be inadequate generally, 
particularly so In the case of thin blades for which an early torsional diver- 
gence was noticed In comparison to MSC/NAS1RAN predictions. It rnay be stated 
that the beam theory used In reference 10 and also the finite element code 
MSC/NAS1RAN, use the geometric nonlinear theory of elasticity In which the 
deflections are treated as large and the strains as small. However, the non- 
linear terms greater than second degree were explicitly discarded In the beam 
theory used In reference 10, while no such simplification Is Implicit In 
MSC/NAS1 RAN. 

Since It has been shown In reference 10 that the addition of second degree 
elastic terms, often referred to as Mil's terms (ref. 11), Is responsible for 
creating an early torsional divergence In the case of thin blades. It Is pro- 
posed to derive the appropriate third degree elastic terms In the coupled flap- 
wise bending, edgewise bending and torsion equations of pretwisted, preconed 
blades. By Introducing the additional, third degree, geometric nonlinear terms 
Into the solution procedure outlined In reference 10, It Is proposed to estab- 
lish that the spurious Instabilities can be removed, and that reliable steady 
state deflections and frequencies can be obtained for complex blade configura- 
tions. Finally, It Is proposed to present and discuss the effects of second- 
and third-degree geometric nonllnearlties on thin and thick blades, and to 
establish the limits of reliability of the present beam theory calculations. 
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EQUATIONS OF MOTION AND METHOD OF SOLUTION 

Figure 1 shows a linearly pretwisted, preconed, and rotating blade of uni- 
form rectangular cross section. The coupled flapwise bending, edgewise bending 
and torsional equations of motion for such a blade with large precone were 
derived in reference 10 including second degree geometric nonlinearities and 
Coriolis effects. By including the additional elastic terms up through third 
degree geometric nonlinearity as given in appendix A, the resulting final equa- 
tions can be shown to be of the following form (a list of notation is given in 
appendix C) : 
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Edgewise bending: 
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Where: 
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It should be noted here that in writing equations (1) to (3), the geo- 
metric pitch angle e p t is replaced by 0 for the sake of convenience in 
writing the equations. This definition of e( = a + yn) in the main body of the 
text should not be confused with the different definition for the same param- 
eter used in appendix A while deriving the curvatures and additional nonlinear 
terms. Furthermore, the additional third degree geometric nonlinear terms of 
elastic origin, represented by H-| , H 2 , and H 3 and in the flatwise, edgewise, 
and torsion equations, are given in appendix A. In the absence of these addi- 
tional third degree geometric nonlinear terms, equations ( 1 ) and ( 2 ) contain 
terms up to 0(c4) in the elastic forces (see ref. 3 for the ordering scheme 
followed for various parameters. For example, w, v", <j> are of the order 0(c); 
w' 4 > is of the order 0 (e 2 ) etc.), and 0 (e 2 ) In the Inertial forces, while 
equation (3), contains elastic force terms up to the order 0 (e 5 ) and inertial 
force terms up to the order 0(e 3 ). It was shown in reference 10 that the addi- 
tion of any further higher order terms, either linear or second degree geomet- 
rically nonlinear, does not alter the quality of the final results. Thus, when 
the present third degree geometric nonlinear elastic terms are added to the 
second degree equations of reference 10 , it is obvious that the order of the 
elastic terms in the bending and torsion equations has increased by at least 
one order of magnitude. For consistency of retaining the terms of appropriate 
magnitudes In the equations, one should, in principle, retain several other 
elastic terms discarded formerly in the second degree equations. However in 
view of the observations made in reference 10 , the terms other than those 
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retained In equations (1) to (3) produce negligible variations on the frequen- 
cies and stability boundaries, and therefore are discarded even In the presence 
of the third degree geometric nonlinear terms. 

Equations (1) to (3) can be written In nondlmenslonal forms by defining 
the following relations: 
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Further, the quantities , H^, and fi 3 in equations (10) to (12) are defined 
by the following: 
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It may be noted here that the last Integral quantities In equations (14) and 
(15) arise due to the consideration of the third degree terms arising from 
foreshortening. Up, In the tension coupling terms (TW)' and (Tv 1 )' In the 
flapwlse and edgewise equations, respectively. The terms that are premulti- 
plied by f 3 0 are the third degree nonlinear terms resulting In the flapwlse 
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and edgewise bending equations respectively from consideration of shearing 
strains In the strain energy variation. Finally, the terms that are premultl- 
plled by f q 2 the two bending equations, and those premultlpl led by f q 3 

In the torsion equation are the nonlinear terms associated with the sectional 
constant, EB^ . The rest of the terms In equations (14) to (16) are the 
nonlinear structural terms of third degree, representing the kinematic pitch 
coupling In the bending and torsion equations (Mil's terms). 

The coupled flapwlse bending, edgewise bending and torsion equations are 
solved by the Galerkln method by assuming the dimensionless deflections In the 
following forms: 
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Where w 0 q , v 0 j , and <j > 0 j are the equilibrium quantities, and Awj , avj, and 
A<t>^ are the perturbation quantities In the generalized coordinates. i|»j and 
0 j are the nonrotating normal mode shape functions for a cantilever beam 
(ref. 12). It may be noted here that sinusoidal mode shape assumed for the 
torsional degree of freedom Is not compatible with the boundary conditions when 
warping Is Included. However, the effect of warping Is not significant for 
large aspect ratio blades considered In this work. Thus, the mode shapes 
assumed here should produce satisfactory results (ref. 10). Proceeding as In 
reference 10, one can apply the Galerkln process for the solution of the non- 
linear steady state equations and the linearized perturbation equations. 
Although It Is not Intended to write all these equations In this work, 
appendix B contains the additional nonlinear terms In the equilibrium equations 
that are over and above those presented In reference 10. These additional non- 
linear terms In the equilibrium equations can be used for writing the corre- 
sponding perturbation equations. 


RESULTS AND DISCUSSION 

The nonlinear steady state equilibrium equations, and the linearized per- 
turbation equations were solved by using computer programs developed In FORTRAN 
language. Integrations are performed on the computer using a Gaussian quadra- 
ture formula. The programs were run on CRAY/XMP computer at NASA Lewis. 

Results were also generated from the finite element code, MSC/NASTRAN, using 
500 CQUAD4 elements, for various cases of rotating, pretwisted, preconed blades 
of various thickness ratios. All these results are presented below. 


Convergence 

The convergence of solutions produced by the Galerkln method with various 
number of nonrotating normal modes In the assumed solution are presented In 
table 1. For the purpose of comparison of the present theoretical results 
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obtained from the third degree equations, Included In this table are the 
results from the second degree equations from reference 10, and those obtained 
from MSC/NASTRAN . In MSC/NASTRAN calculations 250 or 500 CQUAD4 elements were 
used to model the blade. Furthermore, the centrifugal softening effects were 
Incorporated through suitable OMAP/ALTER procedures, solution 64 which accounts 
for differential stiffness effects was used to determine the steady state 
deflections, and subsequently, solution 63 was used to determine the normal 
modes and frequencies. 

From the convergence pattern of the frequency ratios presented In table I, 
It can be seen that a five-mode Galerkln solution produces the lowest five 
coupled mode frequencies that show reasonable agreement with MSC/NASTRAN cal- 
culations, and that the accuracy of higher mode frequencies can be Increased 
by increasing the number of nonrotating normal modes In the assumed solutions. 
Furthermore, a comparison of the results from the third degree equations to 
those from second degree equations Indicates that the accuracy of the third and 
fourth mode coupled frequencies, (corresponding to predominantly first mode 
torsional and third mode flatwise bending frequencies, respectively, that are 
closely coupled for this particular blade configuration) Is improved to a con- 
siderable extent when the present third degree equations are used. 

The convergence pattern of the steady state tip deflections produced by 
the present beam theory is shown In table II. The agreement of present theo- 
retical results to the corresponding ones from MSC/NAS1 RAN is close here also. 

Since a five mode solution (n = 5) Is found to produce the lowest five 
coupled mode frequencies and steady state deflections that are In reasonable 
agreement with the corresponding finite element calculations, further results 
are generated by using a five-mode Galerkin method solution. Such results are 
presented and discussed In what follows. 


Comparison With Experimental Results 

In order to confirm the accuracy of the present theoretical formulation 
including third degree geometric nonlinearities, results pertaining to typical 
untwisted and preconed blade cases of thin rotating blades (corresponding to 
the experimental results reported In ref. 10) are generated. These results 
are presented In table III together with those obtained from the solution of 
second degree geometric nonlinear equations and also those from experiment. 

It should be mentioned here that the elastic modulll used In the theoretical 
calculations are calibrated values (see ref. 10). A mutual comparison of the 
three sets of results presented In table III Indicates that the results pro- 
duced by the present, third degree, beam equations are In closer agreement with 
experimental results than those given by the second degree equations In all 
cases. Considering the torsional mode frequency (mode 2), It can be seen that 
the frequency produced by second degree equations for a 22.5° preconed blade 
with 60° setting angle, and rotating at 3600 rpm. Is 489.9 Hz. Corresponding 
frequency given by the present third degree equations Is 563.8 Hz, while the 
experimental value is 561 Hz. This example Illustrates the nature of the third 
degree geometrically nonlinear elastic terms that are Included In the present 
work. In conclusion, one can state that the results produced by the present 
set of nonlinear equations are In closer agreement with experimental results 
than those produced by the second degree equations. Furthermore, a comparison 
of the results In table III for the particular case of 90° setting angle (blade 
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chord at root section parallel to the axis of rotation) shows that the effect 
of geometric nonl Inearltles for this case of setting angle Is not significant, 
confirming the conclusions drawn In references 9 and 10. 


Vibration and Stability of Pretwisted, Preconed, Rotating Blades 

In order to ascertain the Influence of second and third degree geometric 
nonlinear elastic terms on the coupled frequencies and steady state deflec- 
tions, parametric studies are conducted for rotating, pretwisted, preconed 
blades of various thickness ratio. A typical set of such results are presented 
In table IV(a) for cases of thin blades, and In table IV(b) for cases of thick 
blades. The variation of frequency ratio with the rotational speed Is shown 
In figure 2 for a 15° preconed, thin, blade. 

From an examination of the results presented In table IV(a), It can be 
seen that the spurious Instabilities observed for the case of thin, preconed 
blades, when only the second degree geometric nonllnearltles are Included In 
the equations, are absent when the third degree geometrically nonlinear elastic 
terms are Included to the second degree equations. Furthermore, the quality of 
results Is more accurate when the present third degree equations are used than 
Is the case with those obtained from the second degree equations, even for the 
stable configurations In both sets of equations. 

In order to acquire a further Insight Into the nature of the third degree 
geometric nonlinear terms In coupling the modes, one particular case of rota- 
ting, untwisted, thin blade with 15° precone Is considered. The variation of 
the frequency ratios of the lowest five coupled modes with a variation of the 
rotational parameter, Q/<o-| , Is presented In graphical form In figure 2. The 
curves representing the frequency ratio variations corresponding to predomi- 
nantly flapwlse bending mode frequencies are marked F], F 2 and F 3 ; and those 
corresponding to the first torsion mode and the first edgewise bending mode are 
marked T-] and S-j respectively. A second subscript Is used to Identify the 
degree of nonlinearity retained In the equations that produced the correspond- 
ing results. Thus, the subscript S denotes the second degree equations and 
the subscript t denotes third degree equations. From an examination of the 
results presented In figure 2 , It can be seen that the flapwlse mode frequen- 
cies produced by the second, and the third degree equations are quite close 
(the differences are not apparent In the graphical representation). For this 
particular case of thin, untwisted blade, a clear-cut coupling trend of the 
first torsional and first edgewise modes can be seen In figure 2. When a fur- 
ther comparison of these results Is made to those produced by MSC/NASTRAN 
(shown In figure 2 for discrete values of Q/u-| , by the symbol A), It Is 
observed that the third degree equations produce correct coupling trends, and 
eliminate the torsional Instability that Is shown by the second degree equa- 
tions (refer to the curves marked T-] S and T-|t). The edgewise mode frequency 
predicted by the present third degree equations (S]t) Is closer to MSC/NASTRAN 
results than that predicted by the second degree equations (S-| S ). 

An examination of the results presented In table IV(b) corresponding to 
blades of high thickness ratio Indicates that the second degree equations have 
produced the coupled frequencies that are stable for wide ranges of precones, 
pretwists, and rotational speeds, unlike those In the case of thin blades. 
However, the accuracy of results obtained by using the present third degree 
equations Is seen to be much superior. It Is Interesting to note that while 
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reporting the results In reference 10, it was believed that the Instability 
predicted by MSC/NASTRAN for the case of the thick blade (d/b = 0.2, y = 30°, 
Bp c = 45°, Q/w] = 1.0) was a true static Instability since the second degree 
equations also predicted an Instability at a similar blade rotational speed. 
However, while the present third degree equations are used. It Is found that 
there Is no such Instability at this rotational speed. This observation led 
the authors to re-examine the results of HSC/NAS1RAN, and re-run the code with 
different values for the Inplane stiffening parameters assigned to the Inplane 
nodes of CQUAD4 elements for eliminating the Inplane rotations. It was found 
that the spurious Instability predicted by the solution sequence 64 could be 
eliminated by an appropriate selection of the In-plane stiffening parameter. 

The converged set of results from the present beam theory and HSC/NAS1 RAN are 
thus reported In table IV (b) with an appropriate note being made at the bottom 
of this table. In summary, while MSC/NASTRAN helped In Identifying the limita- 
tions of second degree geometric nonlinear equations, the third degree geomet- 
ric nonlinear equations helped In explaining the spurious Instabilities, 
predicted by MSC/NASTRAN In certain cases. 

From the mutual comparison of the results produced by the second degree 
equations, present third degree equations and the corresponding ones from MSC/ 
NASI RAN, presented in tables I V( a ) and (b). It can be concluded that the spu- 
rious Instabilities observed In the case of thin blades while using the second 
degree equations can be removed by Incorporating the third degree elastic terms 
Into the equations. Furthermore, the quality and reliability of the results 
produced by the present third degree equations are much superior to those 
obtained from second degree equations even In the case of thick blades. For 
configurations that are more complex than those considered In this work, it may 
however be necessary to Include further nonlinear terms, since it appears from 
the present Investigation (refer tables III and IV), that the second degree 
elastic terms have a softening effect (or stiffening effect depending on the 
setting angle and precone) on the torsional modes while the third degree terms 
have a stiffening (or softening) effect on the torsional modes (thus eliminat- 
ing the spurious torsional Instabilities as shown In table IV(a)). While it 
can be contended from the present Investigation that the range of applicability 
of the nonlinear equations Is extended to practical rotor blade operating con- 
ditions by the inclusion of the third degree elastic terms, these are by no 
means, the ultimate form of the equations. Further studies are therefore nec- 
essary for making the beam theory equations applicable for more extreme blade 
configurations. 

Finally, It is believed useful to determine the relative importance of the 
various nonlinear elastic terms of third degree that enter into the equations 
from the consideration of various primary sources In the variational formula- 
tion, such as the normal strain or shearing strains, coupling between tension 
and bending or torsion etc. Such a study will help In Identifying the more 
Important nonlinear terms and In eliminating the relatively less Important 
terms, depending on the type of problem being addressed. To accomplish this 
objective, results are produced by solving the present beam theory equations 
by successively discarding one key effect each time, for a few, typical, pre- 
coned blade configurations. These results are shown in table V, together with 
corresponding results from MSC/NASTRAN. Thus, the results presented under the 
column marked case (1) are obtained from the solution of the third degree non 
linear equations shown by equations (10) to (12). The results presented under 
the column marked case (2) are obtained by Ignoring the third degree nonlinear 
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terms arising from the tension terms (Tw 1 )' and (Tv 1 ) 1 In the two bending equa- 
tions. It may be noted here, by Inspecting equation (4), that these third 
degree elastic terms result from the consideration of foreshortening due to 
bending, up, in the tension expression. Results shown under the column marked 
case (3) are obtained by further neglecting the third degree nonlinear terms 
that arise from the shearing strains. These terms can be Identified In 
equations (A24) and (A25) as those which are associated with GJ. Next, by 
Inspecting equations (A24) to (A26) of appendix A, one can find a group of 
terms In each of these equations premultlplled by the sectional constant, EB^ . 
These terms were addressed In reference 15, and subsequently In reference 16, 
for specialized cases of torsional motion completely uncoupled from bending 
motions of pretwisted blades. It was established that, under the conditions of 
completely uncoupled torsional motion, the nonlinear terms containing the tor- 
sional deformation 4 . and associated with EB-| , are very Important. However, 
their effect In the presence of the kinematic pitch coupling terms Is not com- 
pletely understood. In order to verify the Importance of the nonlinear terms 
associated with the sectional constant EBq In the presence of Mil's terms, 
results are first produced by Ignoring the nonlinear EB-| terms In the two 
bending equations together with those already discarded In case (3). These are 
presented under the column marked case (4). The results obtained by Ignoring 
all the nonlinear terms associated with EB-| In the torsion equation, over and 
above those discarded In the previous case, are presented under the column, 
case (5). Finally, the results obtained by further discarding Mil's terms of 
third degree In the two bending equations but retaining them In the torsion 
equation, are presented under the column marked case ( 6 ). The last column of 
table V shows the results produced by the second degree equations addressed In 
reference 10 . 

From an examination of these results presented In table V, one can observe 
that for the blade configurations addressed In this work, the Influence of 
third degree elastic terms arising from the foreshortening due to bending, up, 
in the bending equations Is of the order of one to two percent on the steady 
state deflections and frequencies. These third degree terms are seen to 
produce a softening effect on the lowest four coupled frequencies. Since the 
corresponding nonlinear terms are not Included In the definition of the ten- 
sion, given In equation (A28) and used for eliminating u', in the torsion 
equation, the exact effect of all these terms Is not known. However, In view 
of the close agreement of theoretical results and NAS1RAN generated results, 
(compare results of MSC/NASTRAN to corresponding ones under cases (1) and (2)), 
It Is believed that the third degree terms arising from tension coupling can 
safely be Ignored In all the equations consistently. Further, from a compar- 
ison of results presented under cases (3) to (5) with the corresponding ones 
from case (2) and MSC/NASTRAN, It Is seen that the third degree nonlinear terms 
arising from shearing strains, and also the nonlinear terms associated with the 
sectional constant, EB-j , produce Insignificant changes on the steady-state 
deflections and on the coupled frequencies. Thus, these nonlinear terms can be 
discarded for blade configurations addressed In this work. Next, the Influence 
of discarding third degree Mil's terms In the bending equations Is seen to be 
quite significant on the steady state deflections and frequencies of pretwisted 
blades. However, the results pertaining to untwisted blades obtained by dis- 
carding third degree Mil's terms In bending equations but retaining In the 
torsion equation are not affected, and remain Identical to those obtained from 
consideration of these terms In all equations. The reason for this would 
become clear by Inspecting the steady state deflections for this particular 
case. It can be seen that for this particular case of untwisted blade with 
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Y = a = 0°, the only significant steady state deflection is w. This gives rise 
to significant equilibrium coordinates and trivial values for v 0 |< and 

<j> 0 |<. Furthermore, since a = y = 0. most of the terms in the equations ( A24 ) 
to ( A26) vanish, and the majority of the still remaining terms of third degree 
will vanish in the perturbation equations due to multiplication of w 0 |< with 
v ok or 4>ok * However, there appears to be a need for caution In discarding 
Mil's terms, since, in the presence of aerodynamic forces (which are not 
addressed in this work), the equilibrium coordinates v^ and may not be 

trivial. In these or similar circumstances, it would be desirable to retain 
third degree Mil's terms in all the equations, whether or not the blade is 
pretwisted . 

After Identifying the necessity and the Importance of the third degree 
elastic terms as discussed above, there remain certain questions concerning the 
rationale of the ordering scheme in discarding the inertial terms. As stated 
in the introduction to this work, and also in appendix A, several nonlinear 
inertial terms were discarded in reference 10 based upon the ordering scheme 
followed therein. For consistency of retaining terms when the third degree 
elastic terms are added to the second degree equations of reference 10, one has 
to consider the influence of the several inertial terms that were discarded 
earlier, before any attempts are made to discard these terms in the presence of 
the third degree elastic terms. In order to answer these questions, at least 
partially, several inertial terms that were discarded previously of order 
0{c 1 * * 4 ), and arising due to rotational effects (which are not Involved with any 
time derivatives), are added to the present torsion equation. On solving the 
resulting coupled equations, it is found that the changes in the frequencies 
were at the fifth significant figure level for a 15° preconed, 30° pretwisted, 
rotating, thin blade (d/b = 0.05, n/u-j = 1.0, a = 0°). In view of this obser- 
vation, it appears rational to discard the large number of higher order 
inertial terms without losing any appreciable accuracy on the final results. 


CONCLUDING REMARKS 

The coupled flapwise bending, edgewise bending and torsion equations of 
dynamic motion of rotating, linearly pretwisted, and large preconed blades of 
symmetric cross section Including third degree elastic, geometrically nonlinear 
terms are derived. These equations are solved by using the Galerkln method and 
a linear perturbation procedure. Parametric results are generated to ascertain 
the necessity and the Influence of third degree elastic terms. Comparisons of 
present theoretical results are made to those produced by the finite element 
code, MSC/NASTRAN, allowing for geometric nonlinear effects; and to those pro- 
duced from an in-house experimental test rig. Close agreement of the present 
theoretical results to the corresponding results from other methods Is observed 
for the parametric range studied. The following specific conclusions have 
emerged In the course of this Investigation: 

(1) For the dynamic analysis of large preconed blades, the second degree 

geometric nonlinear equations are not adequate. The validity of beam theory 

equations can be extended to cover large precones and practical rotational 

speeds by adding the third degree elastic terms involving the kinematic pitch 
coupling to the existing second degree equations. It Is found that third 
degree Mil's terms counteract the Influence of the corresponding second degree 
terms, and thus, the spurious Instabilities that one might observe In the case 
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of thin, rotating blades while using the second degree equations can be sup- 
pressed by considering third degree Mil's terms. It can thus be concluded 
that Inclusion of at least third degree Mil's terms Into the equations Is 
absolutely necessary. 

(2) For the blade configurations addressed In this work, the third degree 
nonlinear elastic terms resulting from the normal strain are found to be Impor- 
tant. The third degree nonlinear terms resulting from the shearing strains 
produce Insignificant changes on the coupled frequencies. Thus, the nonlinear 
terms of third degree that result from the variation of shearing strains may be 
neglected In the equations of motion. Furthermore, the nonlinear terms asso- 
ciated with the sectional constant, EB] , are also found to produce Insignifi- 
cant changes In the steady state deflections and coupled frequencies. In the 
presence of the third degree Mil's terms In the equations, these nonlinear 
terms may also be discarded. 

(3) The present approximation used for defining the tension (only In terms 
of linear parameters In the equation for T) appears to be adequate. Thus, 
defining the tension In all the equations consistently In terms of linear vari- 
ables does not result In any great loss of accuracy In the final results. 

(4) The Influence of higher order Inertial terms, other than those 
retained In the present equations, appears to be Insignificant on the final 
results . 
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APPENDIX A - NONLINEAR CURVATURE EXPRESSIONS AND EQUATIONS OF MOTION 

In order to derive geometric nonlinear equations, It Is necessary to dif- 
ferentiate between the deformed and undeformed configurations, and to Include 
nonlinear terms In the strain displacement relations. If nonlinear terms up to 
third degree are to be consldred In the equations, one has to derive the curva- 
ture expressions having terms up through third degree In the deformation vari- 
ables and also the third degree nonlinear transformation matrix between the 
deformed and undeformed blade-fixed coordinates. In arriving at the deformed 
configuration, one starts with the undeformed configuration and Imposes a 
sequence of rotlonal transformations on the blade-fixed coordinates. Thus, 
deformations u, v, w and 4 > displace the undeformed blade-fixed coordinate 
system x 0 y 0 z 0 to xyz (fig. 3), and rotate xyz to X 3 y 3 Z 3 where X 3 Is 
tangent to the deformed elastic axis. The rotation of the triad xyz to the 
final position may be expressed In terms of the Eulerlan type angles. 

Of the several possible rotational transformation sequences, the one employed 
here Is the flap- lag-pitch rotational transformation sequence. Further details 
on the rotational transformation sequences can be found In reference 13. Thus, 
following an approach Identical to that given In reference 13, and retaining 
terms up through third degree In u, v, w and <$>, the nonlinear curvature 
expressions are derived for a twisted rotor blade undergoing transverse bending 
In two planes (flapwlse along the z-zxls, edgewise along the y-axls), extension 
(u) and torsion (<j>). In what follows, we present the final equations dispens- 
ing with all the Intermediate and lengthy algebra. 


When no assumptions are made on the magnitude of deformation the trans- 
formation matrix between the deformed and undeformed blade fixed coordinates 
Is given by the following equation If nonlinear terms up through third degree 
are retained: 
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The direction cosines given by equations (2) to (10) obey the orthogonality and 
normality conditions up to third degree terms. 

If assumptions corresponding to the case of small deformations I discussed 
In reference 13 are Imposed, and terms up through third degree are retained, 
one obtains the transformation matrix between the deformed and undeformed coor 
dlnate systems as follows: 
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In obtaining equation (11), It Is assumed that the elongations and shears 
are negligible compared to unity, and the squares of derivative of the exten- 
slonal deformation of the elastic axis Is negligible compared to the square of 
bending slopes. The direction cosines given In equation (11) obey the normal- 
ity and orthogonality conditions to third degree terms. It may also be noted 
that the torsion variable, 0 , is the sum of pretwist and elastic twist of the 
blade at any section. Thus, for a linearly pretwisted blade with total pre- 
twist over the length L of the blade, 0 at any distance x from the root 
section Is given by 


0 = a+Yn + <}> = 0pt + 4> 


(A12) 


where a Is the blade setting angle and n = x/L. 


The bending curvatures wy 3 and o Z 3 , and the torsional curvature w X 3 
can be obtained by using either the direction cosines given earlier, or by 
using the projections of the space derivatives of the rotation angles with 
respect to the deformed coordinate S3 along the elastic axis of blade. For 
the case of small deformation level I (ref. 13), the resulting curvature expres- 
sions with terms up through third degree In u,v,w, and (j>, and for flap-lag- 
pitch transformation sequence of rotations are: 
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If r 0 and r-) are the undeformed and deformed position vectors of an arbi- 
trary point on the blade cross section, then the Green's strain tensor, c^j, 
can be obtained from the following equation: 



(A16) 


Proceeding on the same lines as presented In reference 14, one can obtain the 
normal strain y xx and the shearing strains y Xn and y x ^. In terms of the 
curvature expressions, the normal strain y xx , and shearing strains y x ^ and 
Y Xr) can be shown as 
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Considering foreshortening due to bending explicitly as In reference 3, the 
expression for the normal strain to be used In the development of strain energy 
can be shown to be the following: 
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Since It Is of Interest in the current investigation to determine the 
additional third degree elastic terms only, we present here only such terms as 
are obtained from the variation of y xx , y x ^t and y Xn In the strain energy 
variation. Furthermore, we assume the blade to possess symmetry about the two 
principal axes. For this specialized case shown In figure 2, the additional 
nonlinear terms of third degree, over and above those derived In’reference 10, 
are presented below for the flatwise bending (H-|). edgewise bending (H 2 ), and 
torsion (H3) equations, respectively. 
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Additional elastic terms In flapwlse bending equation: 
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It may be noted that the Integral quantities associated with foreshortening, 
up, In equations (A24) and (A25) are the third degree nonlinear terms arising 
from (Tw')‘ and (Tv 1 ) 1 . Furthermore, the underlined terms In equation (A26) 


are not of third degree. 


While the term 




4>‘ 


was considered In the equa- 


tions of reference 10, the term 30' p t <t>' 2 was omitted based upon the ordering 
scheme followed therein. These terms are rewritten here for the sake of avoid- 
ing any confusion. Furthermore, these torsional terms, together with <t>' 3 , 
were discussed by Houbolt and Brooks (ref. 15). Furthermore, It should be 
pointed out here that the effect of higher order terms in the kinetic energy on 
the frequencies and steady state deflections was shown to be of trivial import- 
ance In reference 10. Thus, no attempts are made here to obtain the third 
degree nonlinear expressions pertaining to the variation of the kinetic energy. 


Finally It Is of interest to present the contribution of various terms 
obtained from the variation of the strain energy and pertaining to the exten- 
slonal equation of motion. From the expression of normal strain given In 
equation (A21), one can obtain the following elastic terms belonging to the 
extenslonal equation of motion: 
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From the foregoing expression, one can define the tension T on the blade 
cross section Including third degree nonl Inearl ties as follows by the condi- 
tions of equilibrium of elastic and inertial forces In the extension equation 
(refer to eqs. (95) and (96) of ref. 3): 
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However, while solving the second degree equations In reference 10, only the 
underlined terms were retained In equation (A28) to define the tension, and the 
rest of the second degree geometric nonlinear terms were discarded based upon 
the ordering scheme employed therein. The resulting equation was used to 
eliminate u 1 from the torsion equation. Even with the Inclusion of the lin- 
ear term EA kj( © ' n the tension equation, several nonlinear terms 
result In the torsion equation due to the tension coupling. The effect of all 
such nonlinear terms was found to be negligible on the steady state deflections 
and the coupled frequencies. In view of this observation of reference 10, It 
Is believed that all the nonlinear terms In equation (A28) can be safely dis- 
carded without affecting the final results to any appreciable extent. 
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APPENDIX B - ADDITIONAL NONLINEAR TERMS IN THE EQUILIBRIUM EQUATIONS 

The second degree modal equations resulting from the Galerkln process of 
the coupled bending-torsional equations motion of rotating, pretwisted, pre- 
coned blades were presented In reference 10. The additional nonlinear terms 
that arise In the equilibrium equations due to the consideration of the third 
degree geometric nonl Inearl ties are presented In what follows. These addi- 
tional nonlinear terms are to be added to the corresponding elements of the 
matrix representing the equilibrium equations of reference 10. For this pur- 
pose, the equilibrium equations are defined In the following matrix form. 


Linear 
and 
second 
degree 
nonlinear 
terms 
(ref. 10) 


FI I F2 


El 

i 

1 

E2 

1 

| 

E3 

— - 

i 

-1- 

- — 

4 

— • 

11 

1 

■ 

T2 

1 

T3 


F3 


w 


oj 


iv 


oj 


^oj 
11L J 


= M 


(Bl) 


where FI, F2, .. 13 are the additional nonlinear terms. These elements are 
defined as follows: 


FI -- f,,[w . v „ All,,, „ f d> , 4 ) „ A3T , . 1 * E [v . v „A4T,,. „ * w , w a A5T, .. „ ] 
11 L ok os. 1 jka. M5k oil 1 jka. J s L ok oa. Ijka. ok oa. 1 j ka. J 


+ f 


11 


|y w 


ok v oa. C11 1jka. " 4y V^oa. 031 1 jkfi. * 2y v ok v o9. C41 1 jka. 


* 2y w„,w n C4T 


ok w oa. C4T 1 jka. ' 4y 4> ok 4> oa. c51 1 jka. * 2y V v o! 1 C6T 1 jka. * 2YW ok w oa. C7T 1 jka. 


Y w ok v oa C8T 1 jka. 


]* 


f 1 2 Y v ok v oa. D2T l jka. f f T 0 v ok v oa. D21 1 jka. 


cos 2 0 


fic 


< V 0 kV f W ^ W -o) UF 


ok"o!l ; 


1 jka. 


F2 • f llt V 0 k V 08 A 17 1 jksi - +ok*ol A 2 T 1Jkll f f ll [ YV ok v oa clT 


1 jka. 


( B2) 


" 4y *okV C21 1jkii " Y v ok v o8. C87 1 jka. + 4y 4> ok 4> osi C9T 1jksi 
3a 


2 f 1 2 <t, ok 4, oa D11 ljka. 


F3 =" y f 12 v ok T1 1jk 


(B3) 

(B4) 
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El = 


f n t"ok"oi A71 ijkk - "Wok^ijkkl * f n "ok" 


o9. C81 1 jkS. 


f 4y2 4> ok 4> oe. C9T 1 Jka. " 4y4, ok <t ’o!i C2T ljka, + y w ok w oa. D8T 1 j ka. 


* 2* f 12 4> ok 4> o!l D6T 1jk«. 


(B 5 ) 


E2 ' f ll [ 4y ♦ok*ol C5T 1jkl - 2Y V okV C7T 1jk* * 2y w ok w ol C6T 1Jkl 
- y? v ok w oit D3T 1jk!L * y v okV D4T 1jk*. * 4y2 *ok*ol C3T 1Jkl 
~ v ok v ofi. C41 1 jka. f 2y2 w ok w o8. C47 1 jkS. j ~ f 12 w ok w oi 051 1 jkil 


cos 2 B 


f 10 w ok w o9. D5T 1 jka. + 2 


K < v ok v ok " "okV Uf 1Jkt 


* f ,l t-*ok'*ol A3T 1jkl * v ok"ok A6T 1jk 1 1 f 5 [v ok v ok A8T 1jkk f "okV^ljkt 1 

( B6) 


E3 = y f 1 2 W ok T4 ijk 
T 1 = 0 


T2 - y f -j 3 w ok T8 1jk 


( B 7 ) 
( B8) 
(B9) 


13 71 f 9 [w ok w oa ~ v ok v o!i ] 81 T 1 jks. “ 2f 9 w ok v oa B21 1 jks. 

' 2 f 1 3 *ok T7 1Jk * 3y f 13 w ok v ot D7T 1jk8. _ 3y f 13 4> ok T7 1jk 


(BIO) 


The Integral quantities A1T ! j kSL - - ' D8T i jkS. a PP earing in the e< l uat1ons shown 
above are defined below. 

A ”tJki ■ /|*1 °- 5 *j*k'*i" f °' 25 j s ' n 26 dn 


A2T 


1 jka 


■ / f 4 *j" 9 k*i * 24 , i' e k e i * 2 *j'vi j s1n 20 d ” 
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fl31 ijM = { ♦ 4 *i"v; ♦ 2 »j' e k e i * 2 ^j’ e k e i' | cos 29 <1 '’" 

A4T ijkn ' ( * 4 'j' 4 'k" 4, i * 


x j ^ sin 2 © + cos 2 © | dn 


A5T 1JU ■ J f *lj*j V, 'k'* , i * 4 *ji , k>i" * ,| * , k' , * , i ' J j^2 s1n2 6 * cosZ 9 j dn 

A61 Ukk * ( nj^j'^k"^ * !- 5 'kj'-kk'-ki' * 2 ^"K' * 15 

+ + 0-25 | sin 2© dn 

ATTijkS, = f ^i | 0 - 5 ' ‘•'a, ' f 1,5 * 0,25 s1n 20 dn 

A8T ljki - l * 4 Vk'^" " ♦i*k*i W | |^2 cosZ 6 f s1 " 2 9 j 

A5T t3ki " / - 34, j*k"*i' * ♦j'W 4 4 ’i*k 4, r( 


O O 1 

—z cos^ © *- sin ©> dn 
d d ) 


x |^2 sin 2 0 + cos2 e | dn 


B11 ukt ■ J ( 9 i e j 4 ' i'K' cos 29 d " 

B2T ijki ■ l 9 i 9 j 4, k' 4 'i' s1 " 29 d " 
ciT ijki ■ j T ♦i{*j*k*i" - ♦j*k , *i' j cos 29 dn 


'©. ©„ * 2<jf ! 1 ©. © ' > cos 29 dn 
k 9. j k 9- ( 
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C3T 1JU ■ / Vi 1 Vi cos 26 dn 


C4T uki ' { VjW C0! 20 dn 

C5T uki ‘ ( *i{v"Vi * 2t i' Vi| s1n 26 dn 
C6T uh ’ ( *ijw * *j*k'*i' * °- 5 ♦j'wj s,n 20 d ” 

c,T 'jki • / f 1 - 5 s1n 26 d " 

C8T 1 J kst ■ jf *^*i*k*i’ S '" 26 dn 

"Vt ■ / 


*!©. 0. sin 26 dn 
j k 8. 


‘Vi" d " 


011 1jkt ' l e k e i * 4 *j'°k e i' ” 2 *j 0 k' e i' * 2 *j' 

D2T 'jki ’ f a ’’I j 2>|, j '■‘’k'+i' * * 4 *j" t k'*i f *Y*i*i\ dn 

D31 Uki * jf *l{ 4 *j*k‘*i * s1 " 26 dn 

D4T Ukt * jf 1 *ii 3 *j'*k , *i * 4 *j’kk"'*i * 5 *j*k>i’ * W'| cos 20 d " 
MT 11kt * ( 4, l| 4, i>k' 4 'i' * *W*i"| d " 

D6T 1Jki ’ / ♦l|*j"*k 9 i * 2 *j' 0 M'j d4 

D7i ijk t ■ ( e i \°iK'K' f °j*k"*i * d " 
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D8T ljke. ■ l *1 




cos 20 dn 



T1 uk - X' l 9 W" " 2S W f 9 j"*kJ d ” 

T Sjk ■ /’ *i <e j'" , k' ♦<W"> d " 

T7 Uk ' / VjV dn 



* dn 


The perturbation equations can now be written by Inspecting the equilibrium 
equations presented above and by proceeding as In reference 10. 
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APPENDIX C - NOMENCLATURE 


A 

{B} 

Bl 

b,d 

d/b 

E 

FI ,F2. . .T3 

f 1 » - - * f 1 3 
G 

I.J.k.l 

J 

k A 

k m 

L 

ll . l2» • • - n 3 

m 

n 


0(c) 

P 

R 


r o*n 


s.s 3 

T 


cross-sectional area of blade 
vector 

section constant 

breadth and thickness of blade 

thickness ratio 

Young's modulus 

elements of matrix representing equilibrium equations (see 
appendix B) 

coefficients (see eq. (13)) 
modulus of rigidity 

area moments of Inertia about major and minor principal 
centroldal axes 

dummy Indices 

torsional stiffness constant 

blade cross-sectional polar radius of gyration 

blade cross-sectional mass radius of gyration 

length of blade 

direction cosines 

mass of blade per unit length 

number of nonrotating normal modes for each of the flapwlse 
bending, edgewise bending, and torsional deflections 

small parameter of the order of bending slopes (= w, V or <}>) 

natural radian frequency 

radius of disc 

position vectors of a point on the blade elastic axis, before and 
after deformation, respectively 

coordinate along the elastic axis before and after deformation, 
respectively 

blade tension 
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t 


time 


u,v,w 


uf 

v oj» w oj» < t>oj 

v,w, 4 > 

x 

x 0 yo z o 


displacements of arbitrary point on the elastic axis In x,y,z 
directions, respectively 

foreshortening due to bending (see equation 5) 
steady-state equilibrium quantities 
steady-state dimensionless equilibrium displacements 
running coordinate along x-axIs 

blade-fixed axis system at arbitrary point on elastic axis before 
deformation 


XYZ 


xyz 


y,z 


x 3y3 z 3 


coordinate system with origin at hub (disc) center line which 
rotates with blade such that x-axIs lies along the undeformed 
position of the elastic axis 

blade-fixed axis system which translates with respect to 

x 0 y 0 z o 

centroldal principal axis of beam cross section 

blade-fixed orthogonal axis system In deformed configuration 
obtained by rotating xyz; X 3 axis Is tangent to the deformed 
elastic axis 


a 


Bpc 

Y 


Yxx.Yx n *YxS 


AV j , AWj , A<t>j 


c 


c 1 j 

n 


A ^ 

n.5 


e 


setting angle (collective pitch) 
precone angle 

total pretwist of the blade over Its length 
engineering strain components 
perturbation quantities 

extenslonal component of Green's strain tensor 

Green's strain tensor 

dimensionless length coordinate, x/L 

sectional coordinates along major andmlnor principal axes for a 
given point, respectively 

Eulerlan-type rotation angle In appendix A, eq. (A12), where 
0=oti-7ri^4 > = ©pt * 4> 


6 


geometric pitch angle In equations (1) to (16) and In appendix B 
(e = e p t = a + yh) 
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ej(n) 

Mn.S) 

x n’H 

X 1 

5 

p 

T 

4 > 

'J'j(n) 

Q 

wi 

u x3 » w y3» w z3 

( )' 

(’) 


nonrotating torsional mode shape 
warping function 

derivatives of \ with respect to ri and £, respectively 

frequency parameter, VEI n „/pAL 4 

nondlmenslonal rotational parameter, EI nn /pAI. 4 n 2 

mass density 

dimensionless time, fit 

elastic twist about blade elastic axis 

nonrotating flatwise and edgewise bendlgn mode shapes 

rotor blade angular velocity, rad/s 

exact fundamental mode frequency of straight, untwisted, 
nonrotating beam, 3.51602 X] 

torsional curvature (total rotation about X 3 axis) and 
bending curvatures, respectively 

primes denote differentiation with respect to x or n 

dot over a parameter represents differentiation with respect to 
t or t 
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TABLE I. - CONVERGENCE PATTERN OF FREQUENCY RATIOS (P/Xj) OF PRETWISTED, PRECONED, ROTATING BLADE INCLUDING THIRD 

DEGREE GEOMETRIC NONLINEARITIES AND CORIOLIS EFFECTS 



0.05, 8 


pc 


15 , y = 30 , a = 0 , R = 0, 


L 

d 


200 , - 

0 ), 



Mode 

number 


Results from perturbation solution: Galerkin method 

(Third degree equations) 


Second 
degree 
equations 
(Ref. 10) 
n = 8 

MSC/NASTRAN 
(500 CQUAD4 
elements) 


tn = 1 

n 

B 

B 

n = 5 

n = 6 

B 

n = 8 


6.7837 

4.0123 

3.9880 

3.9820 

3.9813 

3.9802 

3.9804 

3.9805 

3.9699 

3.9923 


65.0281 

40.0635 

20.4535 

20.1407 

20.0567 

20.0226 

20.0116 

20.0073 

19.9823 

20.0707 


74.0808 

63.7502 

60.0351 

57.6998 

57.2708 

57.0335 

56.9803 

56.9702 

51.6155 

57.7069 



89.2954 

85.0776 

63.6660 

63.0034 

62.9864 

62.9480 

62.9221 

60.1018 

64.5911 



216.1802 

137.8281 

86.8464 

86.3469 

86.3378 

86.2964 

86.2841 

85.5334 

86.1357 



438.2073 

214.4627 

214.0674 

124.1565 

121.6039 

121.2614 

120.9963 

120.8147 

121.1038 




351.9987 

275.0052 

214.0269 

204.4178 

200.0460 

199.5745 

199.2196 

199.4268 

8 

— 

__ 

457.2028 

350.2601 

249.5374 

214.1014 

213.9084 

213.9329 

204.6820 

219.5793 


f n represents the number of nonrotating normal modes used in the generalized coordinates. 


TABLE II. - CONVERGENCE PATTERN OF STEADY STATE TIP 
DEFLECTIONS FROM BEAM THEORY INCLUDING THIRD 
DEGREE GEOMETRIC NONLINEARITIES, AND 
COMPARISON WITH MSC/NASTRAN RESULTS 


Method 

Number of 
assumed 
modes or 
CQUAD4 
elements 

Steady state tip deflection 
(centerline deflection) 


w 

t > 

$ 


n = 1 

-0.018444 

0.00182864 

-0.00429201 


n = 2 

-.0547506 

.00661266 

-.00088087 


n = 3 

-.0545014 

.00585071 

-.00086910 

Galerkin 

n = 4 

-.0547382 

.00587672 

-.00081777 


n = 5 

-.0547296 

.00584699 

-.00081642 


n = 6 

-.0547554 

.00584616 

-.00083679 


n = 7 

-.0547519 

.00584231 

-.00083540 


n = 8 

-.0547560 

.00584234 

-.000831394 

MSC / 
NASTRAN 

500 CQUAD4 
Elements 

-.0543090 

.00583720 

-.000902900 



































TABLE III. - COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 


[i 


= 0.05, jj- = 60, y = 0° , L = 152.4 mm (6.0 in), 
p = 2.02728x10”^ kg/mm 3 ] 


E = 222.7 GPa, G = 95.1477 GPa, 


Precone 

angle 

8 pc> 

degree 

Setting 

angle 

O. 

degree 

Rotational 

speed 

a, 

rpm 

Method 

Natural frequency, Hz 

Mode 1 

Mode 2 

Mode 3 

Mode 4 

Mode 5 

0 

30 

3600 

(a) 

126.4 

582.2 

623.9 

1683.4 

1825.8 




(b) 

126.5 

576.2 

625.7 

1689.3 

1813.0 




(c) 


570.0 

630.0 

1688.0 

1807.0 



4800 

(a) 

146.5 

588.7 

647.0 


1837.4 




(b ) 

146.6 

578.4 

651.9 


1815.8 




(c) 

146.0 

582.0 

665.0 

1716.0 

1817.0 



6000 

(a) 

168.6 

597.0 

Mi 

wmm 

1853.3 




(b) 

168.9 

581.5 

ISSI 


1819.9 




(c) 

166.0 

587.0 

B 

1 


22.5 

60 

1200 

mm 

96.5 


593.4 

1655.3 

1796.3 





96.5 

iSH 

593.6 

1656.2 

1796.2 




fail 

97.0 

574.0 

— 


1798.0 



3600 

(a) 

111.5 

489.9 

625.2 

1569.0 

1684.5 




(b) 

114.4 

563.8 

622.9 

1625.7 

1696.7 




(c) 

124.4 

561.0 

626.0 



45 

0 

870 

(a) 

94.8 

570.2 

591.0 

1653.1 





(b) 

94.9 

573.7 

591.3 

1653.9 

1783.6 




(c) 

95.9 

574.0 

588.0 

1655.0 

1799.0 



1506 

(a) 

96.3 

Egg 

593.0 

1655.0 

1675.0 




(b) 

96.5 

WEEm 

593.8 

1657.3 

1696.8 




(c) 

98.6 

Bl 

589.0 

1662.0 

1803.0 

22.5 

90 

2400 

(a) 

SiEI 



Bl 

1815.8 




(b) 

m 89 

in 


ESI 

1815.3 




(c) 

■ ill 

571.4 

602.9 

mm 

1799.0 



4800 

mm 

mm 

573.3 

mm 


1836.1 




ii 

m £9 

573.2 


■ 

1835.7 




mm 

m\m 

562.0 



1826.0 



6000 

mm 

131.8 

572.3 

673.3 


1855.1 





131.7 

572.2 

673.4 

1722.5 

1854.7 





134.5 

553.0 

688.0 

1695.0 

1854.0 


(^Theoretical results excluding third degree geometric nonlinearities (Ref. 10). 
(“'Present theoretical results incuding third degree geometric nonlinearities. 
(“^Experimental results (Reference 10). 






































































TABLE I V ( a) . - COMPARISON OF FREQUENCY RATIOS OF PRECONED, ROTATING BLADES 

OF LOW THICKNESS RATIO 

[a = O’ , R = 0, -jj = 10.] 


d 

F 

Y > 

degree 

B pc, 

degree 

u>l 

Method 

Frequency ratio, p/x^ 

Mode 1 

Mode 2 

Mode 3 

Mode 4 

Mode 5 

0.05 

0 

15 

1.0 

(a) 

(b) 

(c) 

unstable 

5.1562 

5.1781 

23.9231 

23.7926 

42.5942 

44.2514 

64.0009 

63.4078 

90.3709 

88.2218 

2.0 

(a) 

RS 

unstable 

8.2498 

3.3385 

28.5806 

28.3445 

37.0869 

39.0512 

69.5140 

68.0883 

89.9382 

87.6030 

0.05 

30 

15 

0.8 

(a) 

(b) 

(c) 

unstable 

4.6283 

4.6363 





20.4636 

20.5566 

50.6977 

52.6995 

61.8649 

61.7111 

91.1574 

91.3101 

1.0 

(a) 

(b) 

(c) 

unstable 

5.1485 

5.1656 

20.7221 

20.9982 

47.2311 

49.2750 

62.3392 

61.8346 

92.4586 

92.8652 

0.05 

30 

45 

0.3 

(a) 

(b) 

(c) 

3.5517 

3.5652 

3.5762 

19.7350 

19.7493 

19.7675 

54.1394 

57.4509 

57.6886 

60.3105 

63.2751 

65.1148 

84.9439 

85.8371 

85.5516 

0.5 

(a) 

fii 

Unstable 

3.6965 

3.7167 





19.3873 

19.7195 

46.4173 

48.1580 

60.6038 

60.3879 

94.1009 

94.5271 

0.06 

0 

45 

0.3 

(a) 

(b) 

(c) 

3.5438 

3.5541 

3.5665 

22.1084 

22.1433 

22.1613 

50.2454 

52.2279 

52.9289 

61.7771 

61.8699 

61.8449 

72.2320 

74.9998 

74.6194 

0.5 

(a) 

(b) 

(c) 

Unstable 

3.7123 

3.7055 

22.4986 

22.2822 

39.1497 

40.7961 

62.5760 

61.7824 

86.6498 

85.4823 

1.0 

(a) 

(b) 

(c) 

Unstable 

4.8283 

4.8169 

20.6450 

22.2298 

24.8709 

23.0175 

67.6783 

61.1692 

94.9696 

94.4257 


(a)Beam theory including second degree geometric nonlinearities. 
l b |Beam theory including third dgree geometric nonlinearities. 
‘ c ^Results from MSC/NASTRAN. 
































































TABLE I V ( b ) . - COMPARISON OF FREQUENCY RATIOS OF PRECONED, ROTATING BLADES 

OF HIGH THICKNESS RATIO 


[a . 0°, R = 0, = 10.] 


d 

F 

Y > 

degree 

B pc, 

degree 

“i 

Method 


Frequency ratio 

. P l* 1 


Mode 1 

Mode 2 

Mode 3 

Mode 4 

Mode 5 

0.10 

30 

15 

2.0 

■8|§ 

8.2940 

19.6256 

29.0720 

64.0094 

72.7573 






8.2418 

22.1184 

35.5968 

67.5836 

78.0812 





HR 

8.3177 

22.8858 

35.5153 

64.3812 

77.5158 

0.20 

30 

15 

1.0 

(a) 

5.1411 

15.4978 

26.0439 

60.8242 

67.1784 





(b) 

5.1538 

15.5377 

26.3139 

61.3437 

69.0699 





(c) 

5.1619 

15.5424 

26.0745 

60.7458 

70.4579 

0.20 

30 

45 

1.0 

mm 

0.1533 

7.2997 

19.0724 

54.6244 

63.4657 





■Cajfc 

4.8045 

13.1377 

26.5609 

62.3999 

75.6008 





mm 

4.8128 

13.5546 

24.5559 

58.1445 

74.3679 

0.25 

0 

45 

1.0 

mm 

4.6184 

10.6218 

23.2097 

59.4123 

62.8856 






4.8283 

12.2162 

24.8709 

67.6783 

68.7724 





mm 

4.8122 

12.3085 

23.0257 

61.0894 

67.1483 

0.25 

30 

45 

1.0 

m 

4.5877 

10.5075 

23.2236 

56.6582 

63.7352 






4.8175 

11.7024 

25.9483 

61.1655 

71.9986 





mm 

4.8046 

11.8666 

23.9928 

57.2734 

70.3849 


( a )8eam theory including second degree geometric nonlinearities. 
i b 'Beam theory including third dgree geometric nonlinearities. 

Results from MSC/NASTRAN. 

+ Spurious instability was shown by MSC/NASTRAN in reference 10 due to inadequate 
inplane fixing stiffness parameter value. 




















































TABLE V. - EFFECT OF IGNORING CERTAIN GEOMETRIC NONLINEAR TERMS ON THE FREQUENCY RATIOS AND STEADY STATE DEFLECTIONS OF 

ROTATING, PRECONED BLADES 


[a = 0\ R = 0, {j= 10.] 


d 

b 

n_ 

“1 


V, 

degree 

Mode 

number 

or 

tip 

deflection 

MSC/NASTRAN 

Nondimensional frequency parameter, p/x^, ignoring Coriolis effects 

Third degree geometric nonlinear equations 

Second 

degree 

equations 

Case (1) 

Case (2) 

Case (3) 

Case (4) 

Case (5) 

Case (6) 

0.05 

1.0 

30 

15 

Mode 1 
Mode 2 
Mode 3 
Mode 4 
Mode 5 
w 
7 
<t> 

5.1656 

20.9982 

49.2750 

61.8346 

92.8652 

-0.130684 

0.012878 

-0.004148 

5.1485 

20.7221 

47.2311 

62.3392 

92.4586 

-0.132132 

0.013097 

-0.001778 

5.1812 

20.7477 

47.3168 

62.3617 

92.4044 

-0.131552 

0.013022 

-0.001812 

5.1808 

20.7423 

47.3144 

62.3509 

92.4058 

-.0131559 

0.013024 

-0.001808 

5.1807 

20.7420 

47.3102 

62.3515 

92.4084 

-0.131562 

0.013025 

-0.001808 

5.1806 

20.7130 

46.5993 

62.2591 

91.0954 

-0.131558 

0.013023 

-0.001847 

7.0576 

22.5000 

51.5012 

63.1529 

89.9055 

-0.105646 

0.010946 

-0.000827 

Unstable 


0.20 

1.0 

30 

45 

Mode 1 
Mode 2 
Mode 3 
Mode 4 
Mode 5 
w 

"V 

<t> 

4.8128 

13.5546 

24.5559 

58.1445 

74.3679 

-0.357180 

0.035096 

-0.009450 

4.8045 

13.1377 

26.5609 

62.3999 

75.6008 

-0.378955 

0.038377 

-0.006834 

4.9297 

13.2568 

25.5728 

62.4503 

75.3623 

-0.371073 

0.037235 

-0.007093 

4.9265 

13.2127 

26.5587 

62.3605 

75.3583 

-0.371247 

0.037386 

-0.006946 

4.9265 

13.2123 

26.5585 

62.3606 

75.3590 

-0.371248 

0.037386 

-0.006945 

4.9264 

13.2066 

26.5568 

62.3526 

75.2874 

-0.371247 

0.037386 

-0.006942 

5.8022 

13.4365 

27.8104 

61.9055 

74.1130 

-0.321789 

0.037118 

-0.000297 

0.1533 

7.2997 

19.0724 

54.6244 

63.4567 

-0.398513 

0.028370 

-0.069898 

0.10 

2.0 

0 

15 

Mode 1 
Mode 2 
Mode 3 
Mode 4 
Mode 5 
w 
V 
i> 

8.3327 

28.3285 

28.7811 

68.0478 

75.6988 

-0.200225 

8.2498 

27.8202 

28.5806 

69.5140 

77.7080 

-0.205344 

8.4302 

28.0177 

28.7117 

69.6073 

77.4843 

-0.202174 

8.4302 

28.0152 

28.7117 

69.6073 

77.4841 

-0.202174 

8.4302 

28.0152 

28.7117 

69.6073 

77.4841 

-0.202174 

8.4302 

28.0152 

28.7117 

69.6073 

77.4841 

-0.202174 

8.4302 

28.0152 

28.7117 

69.6073 

77.4841 

-0.202174 

8.3781 

21.3722 

28.3088 

68.3325 

70.3606 

-0.203129 


















Full third degree geometric nonlinear equations considered. 

Third degree nonlinear terms arising from (TW')' and (TV 1 )' are ignored. 

Third degree nonlinear terms from (TW)', (TV 1 )' and also those from shearing strains are ignored. 

In addition to those in case (3), third degree nonlinear terms associted with EBx in the bending equations are ignored. 

In addition to those in case (3), third degree nonlinear terms associated with EB^ are ignored in all equations. 

All third degree nonlinear terms are ignored in the bending equations but retained in the torsion equation. Further, all EB^ 
terms are discarded in all equations. 


Case (1) 
Case (2) 
Case (3) 
Case (4) 
Case (5) 
Case (6) 




















































a 




Figure 2. - Comparison of frequency ratios predicted by various methods 
^■Q.05, a -d> =0°, & Pc - 15°, R = oj 

Fj, F^, F 3 are flatwise modes (subscript s denotes solution of second 
degree equations and t denotes solution of third degree equations), 
and Tj and Sj are first torsional and edgewise bending modes, re- 
spectively; A indicates the result produced by MSC/NASTRAN. 
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